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ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ 


Актуальність теми. При регіональних та детальних пошукових роботах в осадових басейнах 
одним з етапів кількісної інтерпретації геофізичних даних є побудова контактної гравітуючої поверх- 
ні. Математично розв'язання таких задач зводять в підсумку до розв'язання оберненої задачі потен- 
ціалу для контактної поверхні. Через використання гармонічних наближень реального розподілу сили 
тяжіння сфера застосування цих задач обмежується областями малої міри з неглибоким заляганням 
джерел. Це звужує клас моделей, для яких можливе адекватне вирішення оберненої задачі. 

В царині контактних задач з часів постановок Нумерова (перша постановка у вигляді неліній- 
ного інтегрального рівняння) та Заморєва (точне обернення оператора рівняння в комплексній області 
за умови, що контактна границя є скінченною ундуляцією прямої на обмеженому інтервалі) користу- 
ються побудовами, які залежать окрім надлишкової густини від середньої глибини до контакту. Заува- 
ження Шванка, яке дозволяло позбутись залежності від глибини П, пройшло повз увагу геофізиків - 
наявна бібліографія свідчить про наближену "лінеаризовану" трактовку контактної задачі у стилі За- 
морєва. Крім того, окреслений Нумеровим клас функцій, на якому розв'язок нелінійного інтегрально- 
го рівняння контакту апроксимується розв'язком лінійного рівняння, не забезпечує однозначного ви- 
рішення задачі. Другий етап розвитку теорії плоскої задачі гравіметрії, зосереджений на досліджен- 
нях в комплексній площині, теж не вирішив усіх проблем, що існують у сфері єдиності розв'язків че- 
рез відсутність (локальних) компактних класів єдиності розв'язків задач з цілком неперервними опе- 
раторами, заданими на нескінченному інтервалі. Пошук обмежень на розв'язки, які б забезпечили їх 
єдиність та стійкість, актуальний і нині. 

Досі теорію визначення контакту двох середовищ з постійною густиною розглядають із вхід- 
ними даними в континуальній області, а на практиці вхідні дані отримують із вимірів певних компо- 
нент поля на істотно коротких профілях. Зростання вимог до точності інтерпретації аномалій спону- 
кає до розвитку відповідного цій практиці теоретичного підгрунтя, а становлення "зрілої комп'ютер- 
ної епохи" в інтерпретації геофізичних полів робить актуальним створення програмного забезпечення 

Зв'язок роботи з науковими планами. Дисертація виконана в Інституті геофізики їм.. С.І. 
Субботіна НАН України у відділі глибинних процесів Землі і гравіметрії в рамках тем , Гранична за- 
дача про відновлення потенціалу за модулем його градієнта та її використання в геодезичній граві- 
метрії 1 в теорії інтерпретації гравітаційних і магнітних аномалій" (Мо держ. реєстр. 01961004773) та 
»Класи єдиності розв'язку обернених задач теорії потенціалу" (Хо держ. реєстр. 01010000236). 

Мета і задачі роботи. Розробка основ "лінеаризованої" постановки оберненої структурної за- 


дачі логарифмічного потенціалу в класі неперервних, разом з першими похідними, функцій для гра- 


З 

ниці поділу однорідних пластів різної густини за умов задання вхідних даних на істотно обмежених 

множинах та розробка дієвих алгоритмів їх чисельного розв'язання в локальній області. Для досяг- 

нення поставленої мети були сформульовані такі основні завдання досліджень: 

1. Визначення інтегральних співвідношень для "лінеаризованої" оберненої задачі структурної 
гравіметрії та класу єдиності її розв'язків за умов збурених похибками вхідних даних у рамках теорії 
наближеного вирішення умовно коректних задач. 

2. Аналіз характерних властивостей "лінеаризованого" інтегрального рівняння для контакту, 
яке дає розв'язок оберненої задачі, узгоджений з даними про похибки вхідних даних та з'ясування не- 
обхідних 1 достатніх умов однозначного розв'язання задачі в аналітичному вигляді за малих збурень. 

3. Реконструкція пропозицій Малкіна, Сенька і Андреєва щодо уточнення наближень Нумеро- 
ва на базі сучасних уявлень про розв'язання некоректних задач. Розробка стійких регуляризуючих 
алгоритмів розв'язання задачі на основі ітераційних процесів Андреєва, Лаврентьєва, Чорного для 
неперервного поля, заданого як на необмежених так і обмежених множинах. 

4. Розробка алгоритмів, позбавлених залежності розв'язку від деякої середньої глибини до по- 
верхні контакту 1 програмного забезпечення для пп. 5-6 та перевірка його якості і надійності на 
тестових даних, вивчення впливу похибок методу і вхідних даних на точність відновлення контакту. 

Наукова новизна отриманих результатів полягає в тому, що: 

е. запропоновано нову постановку задачі для контактної поверхні, яка розширює можливості ін- 
терпретації даних гравіметрії в рамках плоскої теорії потенціалу; 

е вперше отримано ряд нових лінеаризованих рівнянь 1-го роду, які пов'язують задану функцію 
контакту двох необмежених однорідних пластів з вертикальною похідною потенціалу сили тяжін- 
ня. Їх узагальнені аналоги позбавлені параметричної залежності від середньої глибини до контакту 

е. визначено умови коректного вирішення отриманих розв'язків інтегральних рівнянь за різних умов 
задання вхідних даних. Удосконалено способи Малкіна, Сенька та Нумерова-Маловичка для уточ- 
нення лінеаризованих наближень контакту відповідно до постановки задачі; 

е. запропоновано відповідні стійкі швидкозбіжні ітераційні процеси пошуку наближених розв'язків 
рівнянь для контакту на базі теорії узагальнених функцій та еквівалентних перетворень і дослід- 
жено їх збіжність; 

е створено регулярні чисельні алгоритми на основі вищезгаданих ітерацій з типовими обмеженнями 
на параметри середовища 1 поля. 

Достовірність та обгрунтованість розробок підкріплено ретельним і строгим теоретичним 


аналізом математичної моделі і методології вирішення поставленої задачі за допомогою методів фун- 
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кціонального аналізу, теорій потенціалу, узагальнених функцій та розв'язку некоректних задач. При 
викладенні отриманих результатів використано достовірні положення обчислювальної математики, 
методик програмування й інтерпретації геофізичних полів, апробовані результати попередніх дослід- 
жень в цій області. Достовірність результатів обгрунтовано їх співставленням з раніше відомими для 
даного класу задач та апробацією на комплексі тестових прикладів. 

Практичне значення роботи полягає в теоретичному обтрунтуванні лінеаризованого інтегра- 
льного зображення контактної задачі (має загальнометодологічне значення) та в поповненні арсеналу 
технології підбору програмним пакетом регуляризованих алгоритмів для надійного виділення контак- 
тних границь однорідних середовищ за полем, заданим на істотно обмежених інтервалах. 

Особистий внесок автора. Самостійно отримано результати, що описані у всіх розділах, за 
винятком аналізу моделей середовища та ітераційних процесів для обчислення наближень контакту за 


полем, заданим на довгих профілях, дослідження їх властивостей та умов існування розв'язку задачі, 











запозичених з попередніх публікацій |А.В. Чорного, Питання дискретизації рівнянь і побудови регуля- 





ризуючих алгоритмів розвинуто на основі спільних з ним публікацій. Програмна реалізація, тестові 
обчислення і аналіз повністю здійснені автором. 

Апробація результатів дисертації. Положення досліджень, що несуть основне смислове на- 
вантаження та результати чисельного моделювання, доповідались на 28-й сесії Міжнародного семіна- 
ру ім. Д.Г. Успенського "Вопросьт теории и практики геологической интерпретации гравитационньтх, 
магнитньтх и злектрических полей" (м. Київ, 29 січня - 2 лютого 2001 р.), П-й Міжнародній конфере- 
нції для молодих вчених "ГГ еофізичний моніторинг небезпечних геологічних процесів та екологічного 
стану геологічного середовища" (Київ, 8-10 жовтня 2001 р.). 

Публікації. Основні положення дисертації опубліковано в 4 статтях у фахових журналах та 
двох збірках тез (з них одна стаття 1 обидві тези написані автором особисто). 

Структура і обсяг роботи. Дисертація складається із вступу, чотирьох розділів, висновку, пе- 
реліку використаної літератури, бібліографічного коментаря та двох додатків. Загальний обсяг роботи 
200 сторінок, в тому числі 167 найменувань літературних джерел, 30 малюнків та 3 таблиці, обсяг до- 


датків - 26 сторінок. Наукові інтереси автора формувались протягом 1996-2003 р.р. в процесі плідної 











співпраці з науковим керівником |А.В. Чорним, якого автор шанує як вчителя і щиро вдячний йому за 





діяльну участь 1 цінні поради на усіх етапах роботи, що допомогло автору вийти за межі власних при- 
страстей та інтересів. Автор висловлює глибоку вдячність академіку НАН України В./. Старостенку 
за увагу і всебічну підтримку, складає щиру подяку докт. фіз.-мат. наук В.М. Шуману, І.М. Корчагіну 


та канд. фіз.-мат. наук А.М. Тупчієнку за консультації на заключному етапі, а також канд. фіз.-мат. 


б 
наук О.В. Легостаєвій, канд. геол.-мін. наук Т.П. Єгоровій і наук. співр. К.П. Барановій за щире 
ставлення та сприяння у підготовці цієї праці. Окремо виношу подяку сім'ї та любій дружині Фла- 


диславі за розуміння і терпіння. 


ЗМІСТ РОБОТИ 
1. МАТЕМАТИЧНА МОДЕЛЬ ЛІНЕАРИЗОВАНОЇ ЗАДАЧІ ВИЗНАЧЕННЯ КОНТАКТУ 


Як засіб дослідження внутрішньої будови Землі сформульовано математичну модель задачі на 
основі інтегрального подання фундаментального закону взаємодії поля і середовища. Дано детальний 
аналіз наявних математичних моделей оберненої задачі для контактної поверхні. Виявлено, що рівня- 
ння Нумерова і Заморєва, попри зовнішню схожість, математично відтворюють суттєво різні пред- 
метні моделі. Кожне з інтегральних рівнянь, що відповідало тій чи іншій моделі, виявилось однаково залеж- 
ним від двох параметрів: надлишкової густини є і сталої товщини (асимптоти) п шару. Неадекватність та- 
кого опису моделі залишає простір для вдосконалення. За предметну модель задачі обрано найпростішу 
геологічну модель з двох горизонтальних, необмежених, складених породами постійної густини 


с,, і - 1,2 деформованих пластів, що залягають один над одним. Між пластами задано перепад гус- 
тини 0 «с з с, - с. Позначмо через Р, - В. )2р), і-12 необмежені замкнуті області, що містять 


тяжіючі маси й занурені в евклідовий простір Кк). Нехай границі С), цих областей тотожно співпа- 
дають з крівлями першого 1 другого пластів; тоді гравітаційний вплив такої моделі буде еквівалент- 
ним ефекту від області 0 - ), ДЮ, або області Р, з розподілом мас густиною с чи с. Вважаємо, 
що границя 22 - Є, 0 2), області Р кусково-неперервна, її складові ОЇ), - досить гладкі, а сама об- 


ласть досить мала за розмірами, щоб задовольняти умові гармонічної апроксимації потенціалу: 





оо ОМ (х) 9 (хе 
хер | бу дх. 


»0, 


де Ух) -- (о С, р) - збурюючий потенціал сили тяжіння; у - М(х) - внутрішня нормаль до поверх- 


ні рівня (/ (5 со, Ю,) нормального потенціалу / (х) у точці х; мале число є - точність гармонічної 


апроксимації аномалії. Систему координат Ох,, і - 1,3 орієнтовано так, що вісь Ох, спрямована в на- 
прямку тяжіючого шару, а площина Ох,х, щонайліпше апроксимує поверхню Од). Площина х, - 0 
апроксимує поверхню 20), і є її рівнянням на тій підставі, що відомими редукціями завжди можна зве- 
сти виміряні значення сили тяжіння на рівень моря. Нехай складову ОЙ), границі 2) області 0) добре 


апроксимує циліндрична поверхня з твірними, паралельними осі Ох», а функція с), Є є Р розподі- 


6 
лу густин тіл в 0 не залежить від координати 6,. Тоді границею контакту є направляюча 02 «0 (а ) 
2, є Кк? циліндричної поверхні 0/,, яка є простою, без самоперетинів гладкою кривою с(5,). Функції 
контакту відділені від земної поверхні: 0 «ЛП 50! (х) «йртосхаю. 
Для плоскої контактної задачі на основі введеної моделі середовища справедливе рівняння, що 


пов'язує параметри (с, ) середовища й поля 


ге) М 
0, де б 22 "Їм Ге а (1) 


Збудувати модель оберненої задачі для контакту безпосередньо на основі логарифмічного 
потенціалу неможливо через розбіжність інтегралу, який його зображує. Загалом таку задачу зведено 


до розв'язання нелінійного рівняння 





| їн чи - 42 зи(х), и(х)- Р око) 0) (2) 


із значеннями поля, заданими на необмеженій множині. Лінеаризовану постановку контактної задачі 


характеризують два припущення: - коливання контакту малі порівняно з глибиною його залягання; - 


нелінійна частина за своїм вкладом менша від точності спостережень. Клас єдиності розв'язків лі- 
неаризованої задачі характеризує множина Ми (за) (1, р) неперервних за Гельдером разом із своїми 1- 
ми похідними на смузі 0- «од 2 ох «оо й «С(х)ЄЛ,) функцій, які задовольняють умовам 

о(2)« Сх), х є В), 06п «СХ «141, Ок «Ай, 03) 
де сФ(С)єе '- ПО - "коливання" функції ((х), х є ?, й іп Сх), А" «виро (х) - глобальний міні- 
мум і максимум функції 40 (2), - о «жою й є (х) « й". Клас Нумерова МиеУо, р) є декартовим до- 


бутком множини 5 


ос 


(р ) обмежених і (майже всюди) локально інтегрованих функцій, що описують густину 
Се сх, у), (х, у) є Др мас в області Док? та сукупності С"? (в о) Кк «0,,...,; О«а «1 функцій 
увб (х), хе ВО) неперервних (за Гельдером) з їх похідними до К-го порядку включно з показниками а, що 
описують контакти різнорідних шарів зі смуги ІД) -- (х, у): - со «х «оо; 0«С(х)«Є Й, |, Й, - обмежена стала. 
Головну роль в обгрунтуванні апроксимації нелінійного рівняння (2) деяким лінійним оператором на класі 
ми, Р) відіграє така теорема. 

Теорема Нумерова. Для контактних границь, які описуються функціями С (х)х є В"? з 
класу Нумерова Ми ай, І), головну роль в гравітаційному тяжінні підземних мас відіграє перша 


складова у лівій частині рівняння (2), а друга слугує деякою малою поправкою до неї. 
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Висновок 1. У класі функцій Ми (ма)(1, р) нелінійне інтегральне рівняння (2) 1-го роду можна 
замінити з точністю до квадрату ухилення 0! (2, х) ен (2) - сх) лінійним інтегральним рівнянням Фре- 


дгольма 1-го роду 


Де мо), мем) Сх) хє В (осо). 4) 





І і сх) 
я вх) жо) 
Обравши у (4) замість визначеного у фіксованій точці Х значення 0 (х), хє В? будь-яке інше, скажі- 


мо, значення С(у) -п,И «АЙ при у є х, отримаємо нове лінійне рівняння 





Я 2 вух ш -ЛП, У(х)-иїх)- Й. 
іа НО ан с ль Ма)-мд- в (5) 


Рівняння (4), (5) містять повільно спадаючі при І5 -х| -» со ядра Пуассона. Для забезпечення задовіль- 
ної точності розв'язків потрібно задавати функцію и(х) на "майже нескінченній" множині, а виміряне 
поле и(х) отримане на істотно коротких профілях. Для подолання цієї вади є альтернатива: 

-» здійснити високоточну екстраполяцію значень и(х) з короткого на довгий інтервал, після 
чого обчислити наближений розв'язок рівняння; 

-» знайти еквівалентне початковому рівнянню перетворення з швидко спадаючим ядром, яке 
дозволить на коротких інтервалах точно обчислити деяку допоміжну функцію, яка забезпечить визна- 
чення наближення заданого рівняння з гарантованою точністю. 


Поняття істотно короткого профілю введено з таких міркувань. Нехай потрібно визначити в 


рівномірній дискретній мережі точок х, з постійним кроком Дх деякий фрагмент контактної 


поверхні «(х) на інтервалі Їсть ех! за виміряними в інтервалі Їсозез|О Їсть є | з відомою точністю 9, 
ІОх)-чи (2), «б значеннями поля и, (х). Певне наближення фрагменту 2! (х), х є суне | шукають за 


допомогою ітерацій для рівняння (2). Для дискретного алгоритму обчислюють таке значення границі 


інтегрування а, , при якому сумарна похибка апроксимації прямого оператора та похибка чисельного 


т? 
інтегрування не перевищує величини, кратної похибці поля ее п -|ебза)| « 15, І 45. Якщо при цьо- 
му значення поля и, (х) задано на інтервалі Їсозєз так, що Їсо нн М 2 са, ЛК 2| 2 ее то 


інтервал називають майже необмеженим, в іншому разі - істотно обмеженим. 
Для реалізації останнього підходу шляхом перетворень над дійсною частиною інтеграла Шва- 


рца отримано кілька конструкцій з ядрами, що швидко спадають. Наведемо дві з них 














кох п)є ВРУ хе -їк| по ) (6) 


я загіче 
де Д/ - крок чисельного інтегрування. Застосування подібних конструкцій для відновлення контакту 
за умов задання поля на істотно обмежених множинах обгрунтовує така теорема. Нехай 5 (5,1) та 
5 За п)- півсума 1 піврізниця функцій и(ху-П) і и(х, п). 

Теорема 1. Якщо розв'язок й (х) нелінійного інтегрального рівняння (2) для контактної по- 
верхні належить до Мине (І, р), то функція Сх, п) 5 (о, п)-н 5 (бо А), що визначається як сума 


розв'язків лінійних інтегральних рівнянь 1-го роду 


ЕЛ Р Р СО БЕ В ЗМАРННКО сво 
з Їз) ее - ) ази) Р КОЛО 





при і-5, відрізняється від розв'язку С (х) не більше, ніж на квадрат відхилення (оп) а С(х)- В, 
хек? р«н«й. 

Але можна визначити контакт «ра п) на короткому відрізку, виходячи лише з одного рівняння. 

Висновок 2. Якщо розв'язок Сх) нелінійного інтегрального рівняння (2) належить до класу 

ми Є, р), то його лінійне наближення 0 а п) визначається з точністю до квадрату відхилення 

2 (з п) шо (х)-й ,п'«йАА'х є К), як розв'язок лінійного інтегрального рівняння Фредгольма 1-го 


роду: 








У | о РТ (8) 


1 
ла 0 ЦЕ-х)у «(2кХ) 


І Ца л(2- ЗЛ даний збьноді 


4Й 2Й 


Ще один альтернативний шлях для визначення лінійного наближення 0 (зе, п) дає вираз 





Случ 25 ау МЕМ 9) 


де функція 5 "С, п) визначається з першого рівняння системи (7). Ця система забезпечує найточніші 
розрахунки контакту на короткому профілі. Виведено узагальнені аналоги цих виразів. 

Теорема 2. Якщо розв'язок 0 (х) нелінійного інтегрального рівняння (2) для контактної повер- 
хні належить класу Ми" (1, р), то функція Є (х) - 5 (ж) 5 (х), яка визначається як сума роз- 


в'язків лінійних інтегральних рівнянь І-го роду 
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о Її 2) ев Р З 4Є -и(х), 3 їз ємьтесЧає -м), (10) 


відрізняється від розв'язку С (х) не більше, ніж на квадрат відхилення С! (2, х) чі (2)- й (х).Є, хек». 


л(є - х) 








Властивості лінійних операторів системи (10) характеризує така теорема. 
Теорема 3. Якщо функція контакту С! (х) є Ми ФІ, р), то лінійний інтегральний оператор 
из А? Є ) задачі (7) обмежений, неперервний і компактний у просторах С (к0), (вт -1.2.. 
Аналогічний результат справедливий і для операторів (4), (5). Множину допустимих розв язків ліне- 
аризованої контактної задачі характеризують такі обмеження: 
1) контактна границя в п точках, гладка за умовою Та (х)-С (хеп-дх) « оС ) на будь-якому відрізку 
сьо С С,с, З постійним кроком Ах; 


2) у точках профілю х ,, / - 1,2,..М відомі оцінки глибини занурення контакту 0 (х і у» й 


3) апріорі обрано гіпотетичні оцінки 8. (х) подань (х) ро ССх)є сеочде Їсь жа ГЕ; 


4) контактна границя розділяє два однорідні середовища з відомим перепадом густин с є- є - б»; 


5) спостережені в т точках значення поля пс) и(х, )- ди;, і - 1,М містять випадкову похибку за- 


даної інтенсивності Іди, «6; н(х,) є Їсозез 


2. МЕТОДИКА РОЗВ'ЯЗАННЯ ЛІНЕАРИЗОВАНОЇ ЗАДАЧІ ЗА ПОЛЕМ, 
ЗАДАНИМ НА НЕОБМЕЖЕНІЙ МНОЖИНІ 

Викладено методи ітераційного розв'язання лінійних рівнянь для контакту за умов "континуа- 
льного" задання поля та їх узагальнених модифікацій, незалежних від апріорі невідомого параметра Й, 
дано аналіз їх збіжності на основі певних еквівалентних перетворень, наведено оригінальний алго- 
ритм екстраполяції поля при відновленні контакту за традиційною схемою з ядром Пуассона. Довжи- 
ну профілю, необхідну для відновлення з гарантованою точністю аномального контакту, оцінимо так: 
поле на ділянці / - т Дх змоделюємо впливом прямокутника потужністю оС , апріорі більшим за 
ефект досліджуваного контакту; обчисливши із розв'язку прямої задачі, на якій віддалі від цього 
прямокутника його аномальний ефект стане меншим заданої точності є, отримаємо бажану величину І 
Чисельний експеримент доводить, що подовження профілю » на 3/ у кожен бік 1 фіксації кожної 
сублінійної ділянки Дх границі 2-3 точками достатньо для отримання розв'язку із задовільною точні- 
стю в 4,390. 


Вдосконалено класичні способи розв'язання контактної задачі (2). Спосіб Нумерова має низьку 
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точність і забезпечує нульове наближення 0 о) (п), для якого отримано оптимізовані поправки. 
Теорема Андреєва. Якщо розв'язок С (х) нелінійного інтегрального рівняння для кон- 


бі 5 1 щі у 
тактної границі (2) належить до класу Ми! «і, Г), то його лінеаризоване наближення виражаєть- 


ся у вигляді 


З тан Ю й 99 КИ 
(он) з хо о 31)" С 


ла Га а (2 -х) 


у 2 
і від істинного відхиляється не більш, ніж на П (2) при А-» сс. 


|м(2,0)- (0), 1) 





З глибокого аналізу формули (11) отримано наближення контакту, наступні за нульовим: 


(а зх Е Поле з (2 хУЙ и Х и Х -их 
го в) змо) РР уч ху КО» ам) М) мо 
-х)"Їй 


Р З П56/9 -нІЗАЧЕ хх)? (Є 
Сх в) а и(х) Ц 2 Ех ом У 2-х) ну ех 
ЄЗ/о в) з и(х) : ) ізб9бле нІз4ЧЕ - х) 2912-х) -х)'Їй 


Не-ху алу «Є -хУ у нЄ-зУ Цар 


Теорема Малкіна. Корекція кожного з лінеаризованих наближень (12) розв'язку С ех, 1) 








смена: (12) 





рівняння (2) для контактної границі на величину поправки 8 (а) (х, 1) зага Ух, п)ч ЛУ (», 1), де 





ськ Бо З ел) або Аа, (3) 


підвищує точність наступних наближень на порядок порівняно з початковими, тобто до величини 


За 3 
Й (2). Практично достатньо точне наближення поправки Ле (ж, 1) забезпечує величина 


Або вунміоо): синій 


де (Є, ) - деяке "середнє" значення поля на осі Об. 
Теорема Сенька. Уточнення кожного з лінеаризованих наближень (12) розв'язку С Ух, 1) 


рівняння (2) для контактної границі поправкою, яка визначається з лінійного інтегрального рівняння 


першого роду 


І і 2) 
дя (кВ) (2-х) 


генерує наближення С МІРАЛЕ У е ГРЗ ЖАС (х),к - 0,1,2,3..., які ухилятимуться від точного роз- 








- 46 з), Мах) з и(х)- 69(х, п) ета не п 482, а4 
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-3 
в'язку не більш, ніж на величину 25 (р") це тах|с (2). 
Х 
Як узагальнення класичних способів Сенька та Маловичка уточнення контакту запропоновано 
новий ітераційний процес, названий процесом Нумерова-Маловичка: 
жо 2 2 ж 
ПИ Се Б 13 РР й 
Ї гої 5 Р Ї 5 
2 (6-х) юе,(5) я й'к(Є-х) 


ПРИ СБОЖЕНО 
з оз Но) 
сОо)2с5(51), ке0123..те йо, Али. 





АС (хо) - Є, Ох) АСТМІ АС) 


Ах) и(х) - 2, (ж) 46, пабусо, (15) 


Теорема 4. Якщо розв'язок С (х), хе КО уелінійного рівняння (2) для контактної границі на- 
лежить до класу Ми (ад (1, Р), то послідовні наближення 4, (ж, п) є 2, (о А) АС, (ж, п), Я «АХА, 
генеровані ітераційним процесом (15), збігаються при П -» со до граничної функції С (Х), незалежної 
від Й, зі швидкістю геометричної прогресії. 

Розв'язок 0 (х) лінійного інтегрального рівняння (5) відшукуємо у вигляді границі послідовно 


сті га і (з п) , Генерованої ітераційним процесом 





сіна ГОМ ке) с-м обо, 0 б) 


який будемо називати процесом Лаврентьєва- Андреєва. Для спрощення дослідження проблем розв'я- 
зності згаданого процесу ітерацій вдамось до попереднього еквівалентного перетворення виразу (16). 
Лема 1. Якщо послідовність функцій г з (п), генерована ітераційним процесом (16), то ек- 
вівалентна їй послідовність генерується за допомогою послідовних наближень 
с .(слп)а а Їмелиг| сов ОС ЗУ -е"? їі 40. а7) 
Ло 0 ке0 
Теорема 5. Якщо послідовність функцій а б (п) генерується ітераційним процесом (16), 


то вона збігається до граничної функції С (», 1), яка задовольняє інтегральному рівнянню (5) із шви- 


дкістю геометричної прогресії, яка оцінюється нерівністю 


і М 1 и 
Са С5Г) РАДИ - Каз2) Фо кед де ми) з тах|є (о). 








Ця теорема стверджує існування розв'язку 0 (х, 1), а швидкість збіжності спадає з наростанням глиби- 


ни П і номера ітерації. Для узагальненого рівняння (4) існує узагальнений процес ітерацій Лаврентьє- 
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ва- Андреєва, 





сб Го ехмено о) одн лоба 0 аю 


для якого справедливі результати, описані лемою І та теоремою 5, але замість П підставлено значення 
С (х) в деякій фіксованій точці х є ко). 
Для змістовної практичної реалізації в рамках першого постулату альтернативи (значення поля 


1 - є 5 
(0 на майже "нескінченних" інтервалах) запропоновано ефективну схему ек- 


сили тяжіння м(х) хек 
страполяції значень поля за межі профілю, отриману через розклад в ряд Тейлора шуканої функції, 


яка входить до складу інтегралу Пуассона: 





ї 21 21 
0 Т. 1Є----- а шо , 19 
є у ур (ж 2) г - м) (19) 
т 2й 2й 
де мох) - -ЇМе, у удав о| агсів 412. Для відновлення контакту за полем, заданим на 
0 (С -х) 9-х 


майже необмежених множинах, потрібно розв'язати два лінійні інтегральні рівняння Фредгольма 1-го 


роду - (19) та (16) або (18). 


3. ВИЗНАЧЕННЯ КОНТАКТУ ЗА ПОЛЕМ, ЗАДАНИМ НА ІСТОТНО ОБМЕЖЕНИХ МНОЖИНАХ 


Описано методику вирішення лінійних інтегральних рівнянь із швидко спадаючими ядрами, 
призначених для відновлення контакту за полем, заданим на істотно обмежених профілях, проаналі- 
зовано особливості збіжності відповідних ітерацій на введеному класі Нумерова; наведено нову про- 
цедуру прискорення ітерацій. Окреслено множину коректності на основі властивостей відповідних 


лінійних операторів. Так, розв'язок системи (7) шукатимемо за допомогою ітераційного процесу 


Угабьп)ву - зе Їзие -С зе ає чн 85 (ж, 1), по0,со, (20) 


чі ДЕРЕ др), 8 ()ю Зла), 


який надалі називатимемо процесом Лаврентьєва-Чорного. Так само, як і для процесу (16), вдаємось 


до еквівалентного зображення ітерацій. 
Лема 2. Якщо послідовність функцій 15 "бе, І) породжується ітераційним процесом (20), то 


еквівалентна їй послідовність генерується за допомогою послідовних наближень 
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57 (ьк)о ЇМО) и ово - ЗУ(1-( саму). Тло (21) 
У А кед 


Теорема 6. Якщо послідовність функцій 15 М (оп) породжується ітераційним процесом 


(20), то вона збігається до граничної функції 5" (» 1), яка задовольняє інтегральному рівнянню 


1 | тео 5 "(2 пЕ м), 


27 2й 
зі швидкістю геометричної прогресії, яка оцінюється нерівністю 
й М(л) 
СА л)- КУ (» А. Зола де М(л)- тах 5" ( Х, 1). 








з і й з 4 -л 
Прямо можна довести лише збіжність наближень 15 б. (ж и), оскільки нерівність 1- (спам) 


справедлива для будь-яких Л»0 й о»0,а 1-(зват) « Ї виконується лише за обмеження фл» 0.485, 
тому збіжність наближень 15 - (» І) проілюстровано лише чисельно. Співставлення ітераційних про- 
цесів Лаврентьєва- Андреєва (16) і Лаврентьєва-Чорного (20) виявляє, що другий процес збігається 
швидше за перший майже вдвічі, але їх швидкості майже однакові за порядком. Практично процес 
(20) має перевагу над "андреєвським": завдяки швидкому затуханню по горизонталі ядер перетво- 
рень (20) для відновлення з певною точністю наближень 5 (п) потрібно знати значення поля на іс- 


тотно коротших профілях, ніж при використанні інтегралу Пуассона. 


Модуль неперервності оберненого оператора рівняння (20) оцінюється величиною 


Та (х)- кю (2) « тн , де М (л) є тах мох) Усо21) і м(х,2А) - аналітично продовжені значення 


м(х) як гармонічної функції, з рівня у -0 на рівень у - 2. Для відшукання розв'язку узагальненого 


аналогу системи (7) введено узагальнений ітераційний процес Лаврентьєва-Чорного: 


-б ес о3 У РР 
б, ) 5; аа ; ЧУ ає-н8, (х), 


певно у я ЕЛ а ня 2-8 побв, 0 2 


для якого справедливі такі ж твердження, що описані лемою 2 й теоремою 6 для ітерацій (20). Для ос- 











таннього рівняння модуль неперервності має оцінку 
мл") 


б Мв')- 2 птах Р.) (о), 


(х)-є, (а) є 
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де МУ (х) - мажоранта виразу 


руФ(х, СО) - ух) УС 1)" ГОНИЄХ соки УФ(2ає, лє і, и") 


1 
я Око) 2-х) 


за умови (С (х)ея" - тах! .Сх)), а СДР; - аналітично продовжене значення гармонічної функції 
хе 


М (х) зрівня у 20 нарівень у -2//. 


Висновок 3. Кожен з послідовності компактних операторів 


ду; С»х) уні ну Да ех о ар с п -0,1,2,... 





на обмеженій множині функцій м(х) майже проектує на компакт Ми (за) (1, р). 
Звести до єдиного знаменника "швидкісні" характеристики ітерацій (16), (18) та (20), (22) мож- 
на, якщо прискорити збіжність двох перших. Ця процедура поліпшує чисельну стійкість алгоритмів, 


до яких її застосовують. Нехай деяка послідовність наближень 5, (а) є функцією знаменника гео- 


метричної прогресії, 5(4 У 4251 05 (а)- 5, (4) з4",а4 є, власне, знаменником прогресії, на- 





приклад 4 - (лОп-3)) для послідовності 15 З І) ітерацій (20). Для прискорення ітерацій пропо- 


нується така схема 5,(45)- 1, 5, (4,)- 5, (а, Дія, -4; зааИТ 4, 24", «4" . Для ітерацій 


(17) отримаємо "прискорені" наближення у вигляді 


ж 


ські) 2 Реемме | сова хна (о)мо, ау 


-о 


4 (о) 4. (в) з..- 45 (сф), (еф) чо), пе бо. 
Цю схему легко поширити на інші ітераційні процеси. Збіжність наближень (23) описує теорема. 


Теорема 7. Процес послідовних наближень (23) збігається до функції С! (є п) хек, р»д 


м() 


зі швидкістю - и 
НО" -н1 


М(л)- 2 тах|с (х, 1). 





ав ро ага Р ре 


Досліджено стійкість розв'язків поставленої задачі на визначеному класі розв'язків. 


Теорема 8. Якщо за малого є» 0 поля, породжені "паралельними" контактами з класу 


Ми о, р), незначно ухиляються одне від одного в смислі нерівності Ін С)-и, (9| с Є, то самі 


границі незначно відрізняються одна від одної за умовою | р (х)-Є Я (2 сх бЕє,дебеї (аг, п ). 
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4. ОСОБЛИВОСТІ ЧИСЕЛЬНОГО РОЗВ'ЯЗАННЯ КОНТАКТНОЇ ЗАДАЧІ 
НА ОБМЕЖЕНИХ МНОЖИНАХ 


Окреслено загальну стратегію вирішення некоректних задач, здійснено алгебраїзацію пропоно- 
ваних алгоритмів і наведено, у рамках регуляризації за Тихоновим, методику чисельного розв'язання 
лінійного інтегрального рівняння Фредгольма першого роду, до вирішення якого зведено всі дослід- 
жувані процеси. Подано результати чисельного тестування на комп'ютері та порівняльного аналізу 
тестових прикладів, отриманих за різними алгоритмами. Описано програму обчислень на комп'ютері 
наближених розв'язків поставленої задачі на рівномірній мережі вузлів. 


Характеристики прямих операторів АС ; х) задачі 
Абе Г кое МЕ з их), сьожеб, (24) 


дозволяють, в разі точного задання його та правої частини, без проблем будувати таку послідовність 
компактних операторів, яка збігається на певному елементі класу Ми І, р) до точного розв'язку. 


Через нестійкість оберненого оператора потрібно вносити корективи в цю стратегію. 
Скінченновимірну апроксимацію системи (24) здійснено методом проекції. Це дозволяє змен- 
шити розмірність системи, оскільки вхідна матриця симетрична і додатно визначена. В підсумку для 


операторів рівнянь (16) 1 (20) відповідно маємо: 


АЦра)я в, А) а О Кеті, 


АХ (п, х)є Р, (п х)є -- | (аз тезі дети, (25) 


Похибки апроксимації операторів АС з 73 операторами А,, (2 , х) оцінюються, відповідно, як 


«Чі по оаб, |р'бьа, або) 42. 


2-х) з) А 


З цих оцінок знаходять межі інтегрування а, і за формулою Гауса 7-го ступеня точності (12 





є" (з дп ) 





вузлів) здійснюють чисельне інтегрування. По завершенню дискретизації оператора А, г о 6) от- 


римуємо відносно шуканого вектора розв'язку 0 (х, я такі системи лінійних алгебраїчних рівнянь: 


па! 


2п-1 Ах Хо -1 т я т 
на рани 
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2п-2 Ах ау я, б) м и 
з 2. ,.2, со о) Де ДА бо) (26) 


В ході регуляризації за Тихоновим ця система редукується до вирішення рівняння Ейлера 


ак) «УУубо)Кууа, оє Є а) У Кв Ль хо) 028) 


чі ізі 


за умов застосування квадратур Гауса Ї ра (хх лу У у Ї (х уч к( Ж ). Розв'язуючи цю систему для ряду 


значень а: Сб, С- У " , знаходимо шуканий розв'язок 4, КРУ . Оптимізація розв'язків із зна- 
із0 


ходженням а, для модельних прикладів не проводилась. Модельне поле, яке фігурує в правій части- 


ні всіх алгоритмів, отримане внаслідок розв'язання прямої задачі, тобто обчислення поля, генеровано- 


го модельним контактом типу Іванова: 














п 
2 : б 
«Ух її ПРУРНИУН і | СоьС 
Со) й, Д- і їі) 0» Я | (29) 
й (іє со 
Рис 1. Значення отриманого поля, збурені 
же 596-ним шумом, служать вхідними 
з ва даними для вирішення оберненої 
0,2 -е- 5, Сх) пе10 
5 Ох, ), по20 задачі. Під час тестування вивчено 
нен і, С о пн50 . . . 
ол ПРУОРОРІ 5 Г0х, ), по99 чисельну збіжність другого рівнян- 
Р б ня системи (7) (рис. 1) та поведінку 
2- ас зі РАНУ ЧНУ 
оо Рана 5 й - 
и інтегральних ядер (6). Останні нері- 
Р 
9 17" з 7' з ' 3 | 8 в Ж в внозначні за точністю і швидкодією 


Працездатність алгоритмів (16), (13) проілюстровано на модельних прикладах з параметрами: довжи- 


на профілю Іє Їсозе ПЕ І- да 2 км; крок Ах «0.1 км, коливання контакту - 1,6 СС)х -1,4 км; пе- 


репад густини с - 0.3 г/ см); апріорна глибина контакту й - 1,2.км; амплітудні величини а-р- і; 





-(-1Аідх); пе15, Кк, 2-2,1-10, фей екзазеї, Ах « 0133, ДС « 0,066, б - 107. 
Тестування здійснено так: для заданих точних розв'язків 0, (2) іядер К, (ж, 2) розраховують 
аналітично або чисельно праві частини и; (х), потім вносять в них випадкові похибки (нормальний 


розподіл), а точні ядра замінюють наближеними К з (2) тієї ж гладкості з точністю 
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є- Ї ЇК З (х)- К («Р ах .З цими даними вирішено два модельні приклади: 





; Су 1- 2Ї , « 1 
1. Модельне рівняння ГкооеМає з; и(х), срєхос,, точний розв'язок Це) | 5 : г 
ро 
точне ядро к(х, п) -е", при цьому созсз Те |-1.4,1.4|, п -29, Дх - 0.1, се, Ї- І- І, Ц, порядок регу- 
ляризації - нульовий (гладкі функції контакту, правої частини, їх майже однакова ефективна ширина) 
2. Все те саме, але точне ядро К Го п)- дле , що означає внесену похибку оператора є я 0.12 


тобто є 2: 12.59). За такої похибки складно знайти шукан кцію Їх). 
( й ) д укану функцію ((х) 


Ітерації закінчують по досягненні заданої середньоквадратичної нев'язки чи по виконанні за- 
даної кількості кроків ітерації. Похибки отриманих розв'язків не перевищують заданих неточностей 
поля. Порівняння розроблених в роботі 1 відомих класичних алгоритмів Чорної, Нумерова, Малович- 
ка виявило більшу ефективність алгоритмів Лаврентьєва- Андреєва, Лаврентьєва-Чорного порівняно з 
класичними схемами, зумовлену кращим ступенем обумовленості систем алгебраїчних рівнянь. Чи- 
сельні експерименти виявили перевагу алгоритму (16) за точністю над розв'язком (2) Чорної майже 
удвічі, а (20) - майже утричі. З ітераціями на основі рівнянь (8), (9) можна досягти таких же результа- 
тів при вдвічі меншому обсязі обчислень (рис. 2). Дослідження представляють перш за все науковий 
інтерес, але алгоритми можна застосувати у складі автоматизованих систем для моделювання і екс- 


прес-аналізу гравіметричних даних на маловивчених територіях нескладної будови. 
















зоб СО оцоо то сою 
зв - Фо ДО) зо 00 
Рис. 2. я" зт то щу я 0) 
яв з,в ас | - РКО ДР й 
23 3,5- (ї бім 5 
и й -- 39 е х 
о за чи Пе блін 
па із та ги їз є 9.91 
. - с-г ж РТВ Я ом - и 
і, пожити пи їв б" з,2 ЧИ пбтвакь си | 
в си си-ов (ІОВ. Сас ги сис 0 С 
і і і - і і . і і 4 
-7 1 ово С) 1 2 з а 5 6 7 
ві зм 


2-15. ет 
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ОСНОВНІ РЕЗУЛЬТАТИ РОБОТИ 


У дисертації наведено теоретичне узагальнення і нове вирішення оберненої задачі потенціалу 


для контактної границі в лінеаризованій постановці. Головні наукові і практичні досягнення роботи: 


18 
19, Описано математичну модель задачі відновлення контакту на класі Нумерова за вертикаль- 
ними похідними логарифмічного потенціалу у вигляді розв'язків лінеаризованих інтегральних рів- 
нянь Фредгольма 1-го роду. В їх узагальнених модифікаціях усунута залежність розв'язків від апріорі 


невідомого параметра Й -- середньої глибини до контакту. Доведено існування і єдиність розв'язків 


на компактному класі єдиності Мийе(1, р). Модель задачі поширюється на магнітоактивні контакти. 


29. Доведено, що відомі алгоритми, які базуються на виразах з інтегральним ядром Пуассона, 
не забезпечують прийнятної точності відновлення контакту за полем, заданим на коротких профілях, 
які властиві практиці спостережень. Для випадку задання значень поля на істотно короткому профі- 
лі для визначення контактної границі /(х) запропоновано кілька алгоритмів, що базуються на інтегра- 
льних рівняннях із швидкозбіжним ядром типу Шварца. На основі рівнянь з ядрами Пуассона і Шва- 
рца введено два різних типи задач для відновлення контакту (з параметром / та без нього). Вивчено 
структуру класу єдиності зазначених рівнянь. Доведено умовну коректність контактної задачі. 

39. Для розв'язання рівнянь кожного типу задач запропоновано збіжні ітераційні процеси 
Лаврентьєва- Андреєва (16) 1 Лаврентьєва-Чорного (20) та їх узагальнені аналоги (18) 1 (22). Практич- 
ною перевагою ітерацій (20), (22) є швидша збіжність та можливість задання значень контакту на ко- 
роткому профілі. Збіжність усіх ітерацій обгрунтовано за новим способом, що базується на еквіва- 
лентному зображенні підінтегральних виразів в ітераційних процесах через узагальнені функції. 

492, В ході дослідження вищезгаданих ітерацій поліпшено ітераційні способи уточнення нуме- 
рівського наближення контакту, пропоновані свого часу Андреєвим, Малкіним і Сеньком, введено 
новий спосіб Нумерова-Маловичка і обгрунтовано стійкість уточнень. Їх можна використати для ана- 
літичних апроксимацій поля в автоматизованих системах інтерпретації. 

59. Надійне визначення контакту на довгих профілях забезпечене ефективним способом екст- 
раполяції значень поля за межі профілю (на основі розкладу вищих похідних відомого інтеграла Пу- 
ассона в ряд Тейлора), який зведено до розв'язання лінійного рівняння 1-го роду (19). Надійність об- 
числень підвищено ефективним способом прискорення збіжності ітерацій на підставі математичних 
характеристик їх збіжності. 

62. Для чисельного вирішення задачі розроблено регуляризуючі алгоритми на базі концепції 
регуляризації Тихонова. Модернізовано знаходження оптимального параметра регуляризації а, , міні- 
мізовано похибки округлень. Створено пакет програм і перевірено на тестових задачах. Порівняння 
чисельних розв'язків довело перевагу алгоритму (20) над алгоритмом (16). Досягнута точність розв'я- 


зку Є, я 107 при похибках оператора є, е10" і спостережень 0 210" прийнятна для застосувань. 
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льности 04.00.22 - геофизика. - Институт геофизики им. С.И. Субботина НАН Украмньк, Киев, 2004. 

Диссертация посвящена разработке теории и алгоритмов решения задачи восстановления кон- 
тактной границь при условий, что поле представлено значениями производньх логарифмического 
потенциала, заданньжми как на неограниченньх так и на существенно ограниченньх профилях. Ис- 
следовань  основнье математическиє модели обратной задачи логарифмического потенциала для 
контактной поверхности. Отмечено, что параметрическая зависимость решения задачи от плотности 
активного слоя и от средней глубинь или асимптотьт контакта есть следствие рассматриваємьгх моде- 
лей слоистой средь, а не принципиальной характеристикой задачи для неограниченньх тел в отли- 
чие от задачи для тел ограниченньхх. В случає задания поля на длинном профиле предложена зффек- 
тивная схема зкстраполяции значений поля за предельт профиля на основе представления искомой 
функции разложением в ряд Тейлора, а также зкономичная схема ускорения сходимости итерацион- 
ньїх процессов на основе геометрических характеристик их сходимости. Для поля, заданного на огра- 
ниченньх множествах, обоснован ряд предложенньтх способов восстановления контактной границь 
как граничньїжх функций последовательностей решений линейньх интегральньх уравнений первого 
рода с компактньми операторами, имеющими бьттро спадающие интегральнье ядра типа Шварца. 
Моделирование поведения бьїтро спадающих ядер показало, что точность вьтчисления контакта за- 
висит не только от погрешностей поля, но и от длиньг интервала его задания. Установлена корект- 
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ность решения полученньтх уравнений на специально введенном классе Нумерова и разработан регу- 
ляризирующий алгоритм с повьшеной численной устойчивостью. Определен минимальньтій сущест- 
венно ограниченньвй интервал для надежного восстановления контакта, зависящий от соотношения 
длин профиля поля и восстанавливаємого контакта. Осуществлена серия численньх зкспериментов 
для изучения влияния погрешностей на дискретньве аналоги полученньх алгоритмов. Разработанная 
теория без ограничений распространяется на трехмерное пространство и на случай нескольких кон- 
тактньх поверхностей, а также применима и при переменной плотности осадочнь(х слоев. 
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потенціалу на істотно обмежених множинах. -- Рукопис. 


Дисертація на здобуття наукового ступеня кандидата фізико-математичних наук за спеціальні- 
стю 04.00.22 - геофізика. - Інститут геофізики ім. С.І. Субботіна НАН України, Київ, 2004. 

Дисертація присвячена розробці теорії та алгоритмів розв'язання задачі відновлення контакт- 
ної границі за умови, що поле представлено значеннями похідних логарифмічного потенціалу, зада- 
ними як на необмежених, так і на істотно обмежених профілях. В разі поля першого типу запропоно- 
вано ефективні схеми екстраполяції значень поля за межі профілю та прискорення збіжності ітера- 
ційних процесів. Для поля, заданого на обмежених множинах, обгрунтовується ряд запропонованих 
способів відновлення контактної границі як граничних функцій послідовностей розв'язків лінійних 
інтегральних рівнянь першого роду з компактними операторами, що мають інтегральні ядра типу 
Шварца, які швидко спадають. Встановлено коректність розв'язання отриманих рівнянь на спеціаль- 
но введеному класі Нумерова та розроблено регуляризуючий алгоритм з підвищеною чисельною 
стійкістю. 

Ключові слова: логарифмічний потенціал, контактна границя, обернена задача гравіметрії, 
інтегральне рівняння, регуляризуючий алгоритм, ітерація, прискорення збіжності, наближення. 
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